Onset and dynamics of vortex-antivortex pairs in polariton OPO superfluids 
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We study, both theoretically and experimentally, the occurrence of topological defects in polariton 
superfluids in the optical parametric oscillator (OPO) regime. We explain in terms of local supercur- 
rents the deterministic behaviour of both onset and dynamics of vortex-antivortex pairs generated 
by perturbing the system with a pulsed probe. Using a generalised Gross-Pitaevskii equation, in- 
cluding photonic disorder, pumping and decay, we elucidate the reason why topological defects form 
in couples and can be detected by direct visualizations in multi-shot OPO experiments. 

PACS numbers: 42.65.Yj, 47.32. C-, 71.36.+C 



Quantum vortices are topological defects occurring in 
macroscopically coherent systems. Their existence was 
first predicted in superfluids [TJ [2], and later in coher- 
ent waves [3]. Nowadays, quantum vortices have been 
the subject of extensive research across several areas of 
physics and have been observed in type-II superconduc- 
tors, ^He, ultracold atomic gases, nonlinear optics media 
(for a review see, e.g., [H [5]) and very recently micro- 
cavity polaritons [6HT3] . The phase of a quantised vortex 
winds around its core from to 27rm (with m integer) , im- 
plying the vortex carries a quantised angular momentum, 
hm. In contrast with the classical counterpart, quantum 
vortices with the same m are all identical, with a size 
(or healing length) determined by the system nonlinear 
properties. 

Recently, the study of quantized vortices imprinted in 
polariton condensates using pulsed laser fields has at- 
tracted noticeable interest both experimentally [9 and 
theoretically p!4HT7] . providing a diagnostics for super- 
fluid properties of such a non-equilibrium system. In 
particular, resonantly pumped polaritons in the OPO 
regime [18l |19] have been recently shown to exhibit a 
new form of non-equilibrium superfluidity [9", ^20]. Here, 
polaritons continuously injected into the pump state, un- 
dergo coherent stimulated scattering into the signal and 
idler states. An additional pulsed probe can initiate a 
traveling decaying gain, which evolves freely from the 
probe constraints. By using a pulsed Laguerre- Gauss 
(LG) beam, vorticity has been shown to persist not only 
in absence of the rotating drive, but also longer than the 
gain induced by the probe, and therefore to be trans- 
ferred to the OPO signal, demonstrating metastability 
of quantum vortices and persistence of currents [9, 15 . 

However, if the extension of the probe carrying a vor- 
tex with charge m = +1 is smaller than the size of 
the vortex-free OPO signal, continuity of the polariton 



wavefunction requires that necessarily an antivortex with 
charge m = — 1 has to form at the edge of the probe (see 
Fig.[T]). In this Letter, we demonstrate that 'unintended' 
antivortices do appear in the signal at the edge of the 
imprinting vortex probe and explain, both theoretically 
and via experiments, the origin of the deterministic be- 
haviour of the antivortex onset and dynamics. In par- 
ticular, we show where antivortices are more likely to 
appear in terms of the supercurrents of the imprinting 
probe and the ones of the underlying OPO. In addition, 
our study reveals that the onset of vortices in polariton 
superfluids does not require a LG imprinting beam, but 
instead vortex-antivortex (V-AV) pairs can be also gen- 
erated when counter-propagating currents are imposed, 
similarly to what happens in normal (classical) fluids. 

Crucially, via numerical simulations, we elucidate the 
reason why an experimental average over many shots al- 
lows detecting a vortex by direct visualisation in density 
and phase profiles. Recently, it has been suggested by 
stochastic simulations [14 that vortices in non-resonantly 
pumped polariton condensates undergo a random motion 
which will hinder their direct detection, unless they are 
close to being pinned by the stationary disorder poten- 
tial and thus follow a deterministic trajectory [8]. In the 
case considered here of a superfluid generated by the the 
OPO, we can instead explain a deterministic dynamics of 
the V-AV pair in terms of the OPO steady state supercur- 
rents, which determine a unique trajectory for the pair, 
allowing their observation in multi-shot measurements. 

Model The generalised Gross-Pitaevskii equations, 

for coupled cavity and exciton fields ipcxij^^'t) with 
pump and decay [21] model the OPO dynamics {h = 1). 
The exciton-exciton interaction induces a non-linear dy- 




FIG. 1: (Color online) Simulated profile and supercurrents of 
the steady state OPO signal before the arrival of the probe (a) 
and associated interference fringes (b). Location of antivor- 
tices (dots (c)) and vortices (stars (e)) at the arrival (t = ps) 
of a vortex (stars (c)) or an antivortex (dots (e)) probe, for 
1000 realisations of the random relative phase between pump 
and probe, ^rdm- The size of dots in (c) (stars in (e)) is pro- 
portional to the number of times the antivortices (vortices) 
appear in that location. Panel (d) ((f)) shows single shot 
interference fringes relative to the plot in (c) ((e)). Contour- 
level lines in (c) and (e) represent the photonic disorder V{r). 
The white circle represents the edge of the probe. 

namics of lower (LP) and upper (UP) polaritons, the 
eigenstates of the non-interacting Hamiltonian: 

fr _fuJx-ii^x ^i?/2 \ 
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The cavity (exciton) field decays with rate hzc (^x) and 
is replenished by a continuous wave (cw) laser, Fp(r, t) = 
J^fp,cTp{r)e'^^p-''-^p^\ with a top-hat profile [23 . Above 
a pump strength threshold, the system is driven into the 
OPO regime where signal and idler states (with energies 
ujs,i and wavevectors k^^^) get exponentially populated. 
In addition, the OPO is probed by an extra pulsed laser 
Fpi)(v^t). As single shot measurements would give a too 
low signal to noise ratio, an average is performed over 
many pulsed experiments taken always for the same OPO 
conditions. What differs at each probe arrival is the ran- 
dom relative phase ^rdm between pump and probe, 

F(r, t) = Fp{r, t) + Fj,b(r, t)e'*-'"" , (3) 

with ^rdm uniformly distributed between and 27r. 

As already shown in Refs. |9| ^Q, vortices with 
charge m = ±1 can be imprinted in the OPO signal and 



FIG. 2: (Color online) Simulated time evolution of signal af- 
ter the arrival of the vortex probe, averaged over 1000 reali- 
sations of the random phase $rcZm, {il^c{r,t))(^^^^ — spatial 
profile (top) and phase (bottom). Contour- level lines in the 
last panel represent the photonic disorder. 

idler, by adding a LG pulsed probe: 

where the probe momentum kp^ and energy ujpb are res- 
onant with, e.g., the OPO signal state. Here, the pulse 
lasts 2 ps only. The azimuthal angle (p{r) winds from 
to 27T around the vortex core r^^. Finally, to mimic 
the experimental conditions, we include in ([T]) a static 
photonic disorder potential V^(r), with (V^(r)) = and 
{V{r)V{r')) = cr^e-l^-^'l'/^^' {id - 20 /im and aa - 0.1 
meV). 

We solve numerically Eq. ([T]) on a 2D grid by us- 
ing a 5*^-order adaptive-step Runge-Kutta algorithm. 
We first find the steady state stationary conditions for 
OPO emission (/^^ = 0) and plot the OPO signal pro- 
file |V^^(r,t)|e^^^(^'^) by, e.g., filtering in a cone around 
the signal momentum. In addition to the spatial profile, 
|?/^^(r, t)|, we also evaluate the supercurrents V^^(r, t) 
— see Fig. [ija). Note that the presence of the photonic 
disorder does not change qualitatively our results. Its role 
is to break the y —y symmetry left by the pump with 
kp = {kpx^O) and to change accordingly the supercur- 
rents. Further, we simulate the dynamics following the 
arrival of a vortex probe Q at t = ps for 1000 realisa- 
tions of ^rdm and then average (|V^^(r, t)|e*'^^*^'^'^^)<^,^^^. 

Experimental setup The sample studied is a A/2 AlAs 
microcavity with a single GaAs quantum well placed at 
the anti-node of the mirror-confined cavity field, giving 
a Rabi splitting of 4.4 meV - for details on the sam- 
ple see Ref. [22]. Maintaining the sample at 10 K, a cw 
TiiSapphire laser, Fp(r, t), resonantly pumps polaritons 
at 1.5283 eV and kpx = 1.4 /im~^. Above a thresh- 
old, the system enters the steady state OPO regime. For 
a typical pump power of 450 mW, the signal emits at 
1.5268 eV, 1 meV blue-shifted from the LP bare disper- 



3 



Time = 4 ps Time = 35 ps Time = 43 ps Time = 55 ps 




FIG. 3: (Color online) Measured momentum distribution of the vortex probe vs. k — kpb (a). The arrow indicates the signal 
momentum — kpb. Panel (b) ((f)) shows the real space emission of the signal 4 ps after the m = 1 (m = —1) LG probe 
arrival. Panels (c-e) (panels (g-i)) compile, through the interference patterns, the time evolution of an imprinted m — 1 vortex 
(m = — 1 antivortex) and its associated m — —1 antivortex (m = 1 vortex). The straight arrow in (b,f) represents the signal 
current direction in the probe reference frame, while the probe current winds anti-clockwise for m — 1 (clockwise for m — —1). 
The 'unintended' antivortex (vortex) appears in (c) ((g)) at the edge of the probe where the signal and probe currents are 
anti-parallel. V-AV pair dynamics can be followed in both cases for about 30 ps, after which they annihilate. 



sion. We filter the emission in /c-space around the signal 
momentum. A streak camera follows the evolution of the 
signal after the arrival of a 2 ps-long probe pulse from 
a second Ti: Sapphire laser, Fp5(r, t), resonant with the 
OPO signal (typical power 3 /iW), and an average over 
millions of shots is performed. The probe LG profile is 
generated by shining a Gaussian beam through a holo- 
gram with a fork-like dislocation on its fringe pattern. 
Interference images between the signal and an expanded, 
constant phase, region of the signal are obtained in a 
Mach-Zehnder interferometer. 

Results We discuss here the results obtained in the 
numerical simulations and later in the experiments. In 
both cases, we choose OPO conditions such to give a 
vortex- free signal (see Fig. [TJa,b)). Nevertheless, the 
simultaneous presence of pump, signal and idler emit- 
ting at different momenta, as well as the photonic dis- 
order, implies that the OPO steady state is charac- 
terised by currents carrying polaritons from gain- to loss- 
dominated regions. In Fig.jlja), the signal currents have 
a dominant component pointing leftwards and an equilib- 
rium position where all currents point inwards at around 
(-8,-14) /im. 

In single shot simulations of Fig.[TJd,f) (one realisation 
of the phase ^rdm)^ we find that if the probe is positioned 
well inside the OPO signal (e.g., = (0,0) /im), then 
the imprinting of a vortex m = +1 (antivortex m = — 1) 
at t = ps forces the system to generate, at the same 
time, an antivortex m = —1 (vortex m = +1) at the edge 
of the probe. This is a consequence of the continuity of 
the polariton wavefunctions: If the signal OPO phase 
is homogeneous and vortex-free before the arrival of the 
probe, then imposing a topological defect, i.e., a branch 
cut, on the signal phase at the probe core, requires the 



branch cut to terminate where the phase is not imposed 
by the probe any longer and has to continuously connect 
to the freely chosen OPO signal phase, i.e. at the edge 
of the probe. Note that OPO parametric scattering pro- 
cesses constrain the sum of signal and idler phases to the 
phase of the laser pump by 2(j)p = + Thus, at the 
same positions where the V-AV pair appears in the sig- 
nal, an AV-V pair appears in the idler, so that locally the 
phase constraint described above is satisfied. This agrees 
with the experiments in [10 , though there only a single 
V (AV) in the signal (idler) could be detected, because 
the signal size was comparable to the probe one. 

Different relative phases ^rdm cause the antivortex 
(vortex) to appear in different locations around the vor- 
tex (antivortex) probe. However, on 1000 realisations of 
the random phase uniformly distributed between and 
27r, we observe that the antivortices (vortices) are more 
likely to appear on positions where the current of the 
steady state OPO signal before the probe arrival and the 
probe current are opposite. For example, for the m = +1 
(m = —1) probe of Fig. [TJc) (Fig. [TJe)), the current con- 
stantly winds anti-clockwise (clockwise), therefore, com- 
paring with the signal current of Fig. [ija), the two are 
anti-parallel in the bottom right (top left) region on the 
probe edge, region where is very likely that an antivortex 
(vortex) is formed. Note also that the onset of antivor- 
tices (vortices) privileges regions where the steady OPO 
signal has a minimal intensity. Finally, wavefunction con- 
tinuity arguments allows the formation of additional V- 
AV pairs on the probe edge, but are however rare events. 

By averaging the 1000 images obtained at the probe 
arrival (t = ps), e.g., in Fig.jljc), neither the imprinted 
vortex nor the antivortex can be detected (see first panel 
of Fig. [2|: Both phase singularities are washed away by 
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FIG. 4: (Color online) Gaussian probe at rest kpb = (a) and 
moving kpb / (b) shined on the vortex-free OPO signal. 
The measured emission shows that no V-AV pair are created 
if the probe is at rest (c), while a pair appears after about 
15 ps for a moving probe (d). 

averaging the differently positioned branch-cuts. How- 
ever, the steady state signal currents push the V and 
AV, initially positioned in different locations, towards 
the same equilibrium position where all currents point 
inwards. Thus, while at t = ps, on average there is no 
V-AV pair, after ~ 10 ps, both V and AV appear and 
last ^ 75 ps (see Fig. [2|, till they eventually annihilate. 

The theoretical predictions are borne out by the exper- 
imental observations. In order to confirm the role played 
by the relative currents between the probe and the OPO 
signal on the appearance of the 'unintended' antivortex 
(vortex), we inject the vortex (antivortex) probe with 
a finite momentum with respect to that of the signal. 
In Fig. [3](a) we plot the momentum distribution of the 
probe as a function of k — kp^. Thus, in the reference 
frame of the probe, the OPO signal has a definite homo- 
geneous current (straight arrow in Figs. |3];b,f)), while the 
vortex (antivortex) probe has anti-clockwise (clockwise) 
winding constant currents. Figure ^h) ((f)) shows the 
real space emission of the signal 4 ps after the m = 1 
(m = —1) LG probe arrival. Images are taken by sub- 
tracting the steady state OPO. Also, as in Ref. [9 , in 
our experiments, the gain triggered by the probe in the 
signal (and lasting around 25 ps) hinders the observation 
of the underlying signal dynamics, so that in Fig. [3] we 
show the signal evolution after the decay of this extra 
population. According to the previous analysis of Fig. [T] 
we can therefore predict the location of the 'unintended' 
antivortex (vortex) in Fig.jsjc) (Fig.jsjg)), namely where 
the signal and probe currents are anti-parallel. In partic- 
ular, in Fig. [sj^c) the antivortex appears on the opposite 
side of the vortex in Fig. [sj^g). Despite the many-shot 
average, the dynamics of V-AV pairs can be experimen- 
tally followed for about 30 ps (Figs, [sjc-e) and (g-i)), 
thereafter the pair eventually annihilates. 

Finally, we show that it is possible to create a V-AV 
pair with just a Gaussian probe, when there is a difference 
in the signal and probe currents. To this end, we shine 
a Gaussian pulsed beam either at rest with respect to 



the OPO signal, kp5 = kg :^ 0, or moving kp5 7^ k^. 
No pair appears in the first case, while in the second, a 
V-AV pair appears on opposite sides of the probe edge 
— see Fig. [4j Note that circulation is as expected, anti- 
clockwise (clockwise) for the vortex (antivortex) on the 
upper (lower) side of the probe. 

To conclude, the mechanism for V-AV pair formation 
reported here differs from the V-AV binding-unbinding 
associated to the Berezinskii-Kosterlitz-Thouless phase 
transition, recently adopted to interpret the V-AV obser- 
vation in non-resonantly pumped polaritons [11 . In our 
case, the pair onset can be explained in terms of OPO 
and probe relative currents, a simple mechanism which 
does not require resorting to phase fluctuations induced 
by the pump. 
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